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Long-range quantum communication for spin qubits is an important open problem. Here we study 
decoherence of an electron spin qubit that is being transported in a moving quantum dot. We focus 
on spin decoherence due to spin-orbit interaction and a random electric potential. We find that 
at the lowest order, the motion induces longitudinal spin relaxation, with a rate linear in the dot 
velocity. Our calculated spin relaxation time ranges from sub fj,s in GaAs to above ms in Si, making 
this relaxation a significant decoherence channel. Our results also give clear indications on how to 
reduce the decoherence effect of electron motion. 

PACS numbers: 72.25.Rb, 03.67.Lx, 03.67.Hk, 72.25. Dc 



Spin qubits have attracted wide attention over the past 
decade. There have been tremendous progress in the 
experimental study and theoretical investigation of spin 
manipulation and decoherence [l], Q. One of the most 
important advantages of electron spin qubits is that they 
can be coupled strongly via the exchange interaction |3| , 
which allows fast two-spin gates. However, exchange in- 
teraction is short-ranged, and long-range quantum com- 
munication remains a significant open problem in the 
scale-up of spin qubit architectures. 

Various ideas have been proposed to move spin infor- 
mation on chip, such as via spin-photon coupling 0, 
spin bus @ , and directly moving the electrons themselves 
[7H10|. The last approach is attractive because of its 
conceptual simplicity and its similarity to the conven- 
tional charge-coupled devices. In recent years, several 
experimental groups have shown how a surface acoustic 
wave (SAW) can transport an electron over a distance of 
several microns between quantum dots 11- lB|. Clearly, 



transferring spin information by directly moving its car- 
rier is an intriguing and promising approach, and de- 
serves further in-depth analysis. Here we focus on the 
aspect of spin decoherence due to electron motion. 

The main decoherence channel for a confined electron 
spin in a finite field is the hyperfine interaction induced 
pure dephasing [l6|, 17 1. Spin relaxation due to spin- 
orbit interaction is much slower [3]. On the other hand, 
spin relaxation of free electrons and holes in semiconduc- 
tors is dominated by spin-orbit interaction , while 
the effect of hyperfine interaction is strongly suppressed 
by motional narrowing (20| . For a moving electron spin 
qubit with controlled motion, an intriguing question is 
thus when decoherence due to spin-orbit interaction be- 
comes dominant. 

In this Letter we study spin decoherence of a moving 
but confined electron due to static disorders in a semicon- 
ductor heterostructure. For example, in a modulation- 
doped GaAs/Al^Gai-^As structure, the ionized dopant 
nuclei produce a random electric potential at the GaAs 
interface where the quantum dot (QD) is located. If the 
QD is moved (by programming the top gate potentials) 



along the interface, the confined electron would experi- 
ence random kicks from this potential. Through spin- 
orbit interaction, the electron spin can sense this ran- 
dom potential as well, which leads to spin decoherence. 
Here we construct a theoretical description of this deco- 
herence mechanism, and calculate the corresponding spin 
decoherence rate. We find that this is a longitudinal re- 
laxation channel at the lowest order. Its rate can be as 
fast as sub-/xs in GaAs or above ms in Si, making it a 
significant decoherence channel. 

Electron Hamiltonian in the Moving Frame — The 
model system we consider is a single electron in a gate- 
defined QD from a two-dimensional electron gas (2DEG). 
In general the growth-direction confinement is much 
stronger, so that the QD can only be moved perpendicu- 
lar to the growth direction by programming the top gate 
potentials. The electron is always confined while the QD 
is moved. Indeed, we assume the QD motion is suffi- 
ciently slow (adiabatic) so that the electron remains in 
the ground orbital state. 
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FIG. 1: A schematic of a spin qubit in a moving QD. Panel (a) 
gives the topview of the structure and the coordinate system 
(xyz) defined in the laboratory frame, with x and y along the 
[110] and [110] directions. Panel (b) gives the sideview and 
the total effective magnetic field. 

As shown in Fig.Q] we consider a uniform linear motion 
of a QD with a constant velocity Vo (the QD potential 
minimum is at ro (t) = v^t = [xq (t) ,yo (t)]), for which 
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it is convenient to use a moving reference frame (coor- 
dinates x'y'z', while the coordinates in the laboratory 
frame are xyz, with z and z' along the growth direction) 
of the same velocity. The system Hamiltonian in the 
moving frame is given by (2lj 



tain the effective spin Hamiltonian [2l| 
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Here r = [x,y] (r' = [x',y']) and tc (it') are the elec- 
tron two-dimensional position and momentum operators 
in the laboratory (moving) reference frame. Operators 
in different frames are related: r = r' + ro (t) with 
ro (t) = Vot. The sub-indices d, Z, and SO refer to 
the dot, Zeeman interaction, and the spin-orbit inter- 
action, and the random electrical potential is captured 
by 5V (r). In the moving frame, the random poten- 
tial is time-dependent due to the QD motion, SV — 
SV (ro (t) + r'), so that the static disorder is now a 
charge noise. In Hd, n' = — iKV' + (e/c)A(r') —pL—Pa 
with p L (t) = -eB 0z /c[-y (t) ,x (t) ,0], p = m*v 
and [^jTTj] = —ifi^£ijkBok, in which Bq is a constant 
and uniform applied magnetic field. U (r') = |m*Wjr' 2 
is the time-independent harmonic confinement potential 
in the moving frame. In Hz, B = Bq + AB is the total 
magnetic field in the moving frame, in which AB is the 
orbital motion induced effective magnetic field 

AB = ^ m l ° sin <p v ,/3 + cos </>„, 0) , (5) 
3M-B 

where <p v is the angle between the dot velocity and 
the [110] crystal axis. The Zeeman frequency is ljz = 
g^sB/h and the spin quantization direction is n = B/B, 
which is generally different from the direction of the ap- 
plied magnetic field Bo- In Hso, P± = (fi ± a) are the 
coupling constants, where a and B are the Rashba and 
Dresselhaus coupling constants. 

Effective Spin Hamiltonian — With the QD motion adi- 
abatic, and the electron orbital degree of freedom not 
excited, we can focus on the electron spin dynamics 
by decoupling the spin space (with the ground orbital 
state) from the rest of the Hilbcrt space. Specifi- 
cally, we perform a Schricffcr- Wolff transformation H = 
cxp(S)H exp(— S) to remove the spin-orbit coupling 
(SOC) in the leading order [3, [22m3|. This requires 
[Hd + Hz , S] = Hso , which means S ~ icr ■ £ + 
hx> z (n x L^ 1 ^) • er, assuming that m* (/3 2 + a 2 ) <C 
huz *C hujd- Here L is the Liouville super-operator, 
hdA = [Hd,A], and £ = (y'/X-, x'/X + , 0) is a vector 
in the 2DEG plane, where A± = Tij (m*j3±) are the spin- 
orbit lengths. After simplifying the commutators, we ob- 



H eff 
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n{t) 



= ^ gl x B (B + SB(t))- a 
= 2B x $1 (i) 
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This effective Hamiltonian has two important features. 
First, the spatially random electric field e(r) is now a 
temporally random magnetic field for the electron spin, 
SB(t). This transformation is through the dot motion 
ro (t) [e(r) — > e(t)} and the spin-orbit interaction [e(t) — > 
5B(t)}. Second, Equations (JT]) shows that in the first 
order in SOC, there can be only transverse fluctuations 
in the effective magnetic field since SB(t) ■ B = [See 
Fig. Hlb)]. D ue to this transverse nature, there is no pure 
dephasing in the first order of SOC. 

Noise Correlation — To calculate the spin decay due to 
5B(t), we need to first calculate its correlation functions 
(SBi^SBj (t)). As illustrated in Fig. [TJ when the QD 
moves, the field 5B (t) fluctuates as a result of the SOC 
and the momentum scattering by the random potential. 
Thus the temporal correlation of the random field SB (t) 
is directly related to the spatial correlation of the random 
electric field in the laboratory frame (ej (ri) Sj (r 2 )). 

We assume that the random electric field (due to mod- 
ulation doping, or other disorder) is isotropic, 



(ei (n)sj (r 2 )) = S., 



-Ar/l e 



where 5ij is the Kronecker delta function (i, j = x' or y'), 
a £ is the standard deviation of the electric field, Ar = 
— 7*2], and l £ is the correlation length of the random 
field. Here we assume that the cut-off function takes the 
exponential form exp(— Ar/l £ ), but it can also take on 
other forms f c (Ar/l £ ) [2l| . To simplify the expression 
of the spin decay rate, we rotate to a new coordinate 
system XYZ, where Z is along the total magnetic field 
direction h, so that SB is in the XY plane. The XY 
axes are further chosen such that (SBx(Q)SBy (t)) = 0. 
The magnetic field correlation can then be calculated as 
21 



(SB x {0)SB x (t)) = 



2eBa £ 



(9) 



where t c = l s /vo is the correlation time. The effective 
SOC length A + is given by 



1 
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by substituting 



(SBySBYit)) is obtained from Eq. 
A_ for A+. 

Spin Relaxation — Now we study decoherence of the 
electron spin qubit due to Hamiltonian ([5]). The noise 
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correlation time r c is generally much shorter than the 
qubit decay time (the inset of Fig. [2] shows values of r c ). 
In this regime, the dynamics of the spin qubit are gov- 
erned by the Bloch equations [25|. With purely trans- 
verse fluctuations, the longitudinal and transverse relax- 
ation rates, 1/Ti and 1/T 2 , are [3, [H, HH 



1 



2 

2h 2 



J xx( U} z) + Jyy^z), 



(10) 
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(SBi^SBjit)) cos (wt)dt. 



Using Eq. © and its Jyy correspondent, we obtain 



1 

Fso 



2u 2 z t c 



so\ 



J Z'c 



1 



0), (11) 

2a(3 sin 2 cos 20] .(12) 



Here and ^ are the polar and azimuthal angles of B in 
the x'y'z' coordinates. 

Before delving into the numerics we first discuss some 
qualitative features of the spin relaxation rate here. 
Firstly, 1/Ti oc 1/w^- This strong dependence on the 
QD confinement means that this spin relaxation channel 
can be suppressed by having strong confinement for the 
QD. Secondly, 1/Ti oc a 2 . The origin of the static disor- 
der would determine the magnitude here. For example, 
in a modulation doped GaAs structure, 5V ~ 20 mV [26( 
and l e ~ OT^m j2f|, so that a £ = SV/l s ~ 200 kV/m. 
On the other hand, for an undoped top-gate structure 
in Si 27 1, there could be disorder from defects in the 



barrier, though its characteristic length and strength are 
unknown (most probably much smaller than in the modu- 
lation doped structures) . Our numerical estimates below 
use parameters from the modulation doped structures. 

The SOC-dependence of 1/Ti is contained in Fso m 
terms of a and f3, the Rashba and Dresselhaus SOC 
strength. These parameters are materials- and device- 
specific. In Si /? = 0, while in GaAs floaAs = 300 m/s 
is fixed [2l[. In both materials a is dependent on the 
particular quantum well structure and doping. 

The dependence on the direction of magnetic field B 
by 1/Ti is also contained in Tso, in terms of the polar 
and azimuthal angles 8 and (f>. For example, for a per- 



pendicular field (B 
in-plane field, Fso 



[001]), F so = 2 ([3 2 + a 2 ). For an 
[3 2 + a 2 + 2a[3cos2(b. Thus, the 
decay rate in a perpendicular field is always larger than 
if the field is in-plane (1/Ti) perp > (1/Ti) in _ piarie . _ 

In the case of an in-plane magnetic field, the spin re- 
laxation rate 1/Ti has a sinusoidal dependence on the 
azimuthal angle <fi of the _B-field. The minimum rate 
is 1/Ti = 2 [ea e [fi - a) / (fejg)] 3 w|r c / (l + oj 2 t 2 ) (as- 
suming a/3 > 0), when the B-field is along the y' axis 
(<f) = ir/2). Thus, in the special case when a = f3 
and 4> = it/2 (or a = -(3 and 4> = 0), 1/Ti = 0! In 
other words, since AB is along y' (x') axis when a — f3 



(a = —0) (c.f. Eq. ([5|)), spin relaxation due to QD mo- 
tion vanishes if the applied magnetic field Bq is along 
y' for a = f3 (or along the x' axis for a = ~j3). Such 
special cases (a = ±/3) have been discussed previously 
in the context of spin relaxation due to phonon emis- 
sion [l8|j [1^] . Note that Hamiltonian H in |T]) conserves 
the spin component oy^') for a = /3 (a = —j3) and 
B || y' (x 1 ). This spin conservation results in Ti being 
infinite to all orders in the spin-orbit Hamiltonian (j4|). 
Meanwhile, decoherence rate I/T2 reduces to the next 
order contribution of in the form of pure dephasing. 
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FIG. 2: The motion induced effective magnetic field AB as 
a function of moving velocity for GaAs (solid line) QD with 
13 — 300 m/s and Si (dashed line) QD with a = 5 m/s. The 
inset gives the bath correlation time r c . 

The dependence on the magnitude of B and velocity 
v by 1/Ti is contained in uj 2 z t c / (l + w|r 2 ) of Eq. (JTJJ) . 
Here we first estimate the magnitude of the motion in- 
duced field AB. In Fig. [2 we plot the motion induced 
magnetic field AB as a function of velocity |t>o| for GaAs 
and Si QDs. The magnitude of AB in GaAs is two or- 
ders of magnitude larger than in Si, but still negligible 
if a strong magnetic field (order of Tesla) is applied. We 
discuss the high and low field cases separately as follows. 

High field and slow moving limit — For a strong applied 
magnetic field (B > 1 T) and a slow moving QD (1 nm/ns 
< «o < lOOnm/ns), AB can be neglected, and the con- 
dition u)zt c 3> 1 is satisfied. In this limit, the loz (or B) 
dependence cancels out in the Eq. pT|) . 1/Ti depends 
linearly on the speed Vq of the QD, and is independent 
of the direction of the motion (since AB is neglected) . 

We carry out numerical calculations on two representa- 
tive QD structures, one in GaAs/Ali-^Ga^As, the other 
in Si/SiGe. In both cases, the dot confinement energies 
are set at fmJd = 1 meV and 3 meV, and the applied 
magnetic field is B = 1 T. For the GaAs QD, we use the 
bulk g-factor g = —0.44, and the electron effective mass 
m* = 0.067mo, where mo is the free electron rest mass. 
For the Si QD, g = 2, m* = 0.19mo, and the Rashba 
SOC strength is chosen to be asi = 5 m/s [2§-31|. Fig. [3] 
shows the spin relaxation rate 1/Ti as a function of the 
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FIG. 3: Spin relaxation rate 1/Ti as a function of the QD 
velocity for GaAs QDs with j3 = 300 m/s and Si QDs with 
a = 5 m/s (with an in-plane magnetic field). Here ojd and SV 
are in units of meV and mV. 



QD speed in an in-plane _B-ficld, when 

-plane ^ e 



i 



(/3 2 



a 



2a/3cos20) . 

(13) 

For a moving GaAs QD, we find that Ti ranges from 
/xs to 10 ms. For a Si QD, T\ > 1 ms because of the 
weaker spin-orbit coupling. In terms of the moving dis- 
tance (voTi) in _ plane , spin coherence is lost in as short as 
3 /jm in GaAs and 0.1 m in Si. Clearly, while the spin 
relaxation times here are still much longer than those in 
a 2DEG, the QD motion does present a serious threat to 
the coherence of the spin qubit, especially in modulation 
doped GaAs heterostructures. 

Low field and fast moving limit — If the magnetic field 
is low, and/or the QD motion is fast (but still adiabatic), 
so that LUzT c <C 1, the spin relaxation rate is 



1 



F so {9,(l>W 2 



(14) 



If the applied field Bq is much larger than the motion 
induced field AB, so that g^BAB/h <C loz <C t" 1 , 
we obtain 1/Ti V^Oj indicating motional narrow- 
ing. Whereas, if B = 0, only the motion induced field 
AJB contributes to the spin splitting, with 9 = ir/2 and 
/?_ tan<^> = j3 + cot <j> v . Noww|r c = (2m* /H) 2 vqI e F v (<f) v ), 
where F v (<f) v ) = (f3 2 + a 2 + 2a(3 cos 2<f) v ) . Interestingly, 
this (/^-dependence is completely canceled by that in 
Fso , so that 

1 r ,i2 



= 2v l e 



2ea e m* (fi 2 - a 2 ) / (huj d )' 



(15) 



In other words, when no magnetic field is applied, 1/Ti 
depends linearly on the speed vq of the QD motion and 
is independent of the direction <j) v of the motion. 



As an example we consider an SAW-confincd electron 
in GaAs. Here the QD moves fast, at the speed of 
sound v = vsaw = 3000 m/s, so that the low-field 
limit ijOz t c < 1 is satisfied even for B ~ 1 T. The 
electron should still remain in the ground orbital state, 
however, while the motion-induced magnetic field is now 
AB ~ 0.02 T. The spin relaxation rate should have a 
weak dependence on the direction of the motion when 
Bq and AB are comparable. The confinement potential 
for an electron in an SAW is dependent on the driving in- 
tensity PsAWi with ujd ~ 1 meV [2l|. Using parameters 
for modulation doped samples, l £ ~ 0.1/xm and a e ~ 200 
kV/m, we estimate that 1/Ti ~ 10 8 Hz in an strong 
in-plane magnetic field (B ^> AB, f3 = 300 m/s and 
a = 0) and l/T t - 10 5 Hz when B = 0. These rates 
can be reduced by having a larger uid (with higher SAW 
driving intensity and / or higher frequency) , and most im- 
portantly, using a less disordered sample with smaller a e . 

Wc note that the energy exchange involved in the spin 
relaxation studied here is between the electron spin and 
the external agent that drives the QD motion. In other 
words, we are dealing with an open system and energy 
is not conserved. It is also interesting to note that if we 
reverse the direction of the QD velocity, the motion in- 
duced field AB reverses its sign, while the magnetic noise 
SB does not (Bq ^ 0), so that its effect is irreversible. 
However, if there is no applied magnetic field Bq = 0, 
the motion induced field AB and the magnetic noise SB 
are both reversed if the direction of the QD velocity is 
reversed. The energy exchange with the external agent 
can be exactly canceled out in this case. 

In conclusion, we have studied spin relaxation of an 
electron confined in a moving QD. The spin relaxation 
mechanism we studied originates from momentum scat- 
tering and spin-orbit interaction. At the lowest order, it 
is a longitudinal relaxation channel, so that = 2T\. 
The relaxation rate is inversely proportional to the 4th 
power of the confinement energy, so that the spin de- 
coherence is the fastest for large quantum dots. Fur- 
thermore, under the condition of high-field slow-motion 
or very-low-field fast motion, the decoherence rate in- 
creases linearly with the speed of the QD. Quantitatively, 
wc find that in modulation-doped GaAs heterostructures 
this can be an important spin decoherence channel, where 
the spin relaxation time can be as short as sub-^s and as 
long as ms, depending on the QD confinement strength 
and the magnitude of the random potential. For mod- 
ulation doped Si/SiGe QDs the spin relaxation rate is 
generally several orders of magnitude slower. However, 
compared to known spin decoherence channels in Si, this 
relaxation can still be quite significant. 

We thank support by US ARO (W911NF0910393), 
DARPA QuEST through AFOSR, and NSF PIF (PHY- 
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Supplementary electronic material for: 
Spin Qubit Relaxation in a Moving Quantum Dot 

In this supplementary we give an expanded discussion on how some of the results in the main text are obtained. 
We first derive the Hamiltonian in the moving frame as used in the main text. We then show how to obtain the 
effective spin Hamiltonian for the moving electron. We also derive the correlation functions of the effective magnetic 
fluctuations. Finally, we discuss the case of a surface acoustic wave (SAW) trapped electron spin. 



TRANSFORMING TO THE MOVING FRAME 



The key to the study of decoherence of a moving electron is to disentangle the spin dynamics from the orbital 
dynamics. In the first step to achieve this separation, we transform into the moving reference frame in which the 
electron is at rest in the absence of environmental randomness, and the QD confinement potential becomes time- 
independent. The specific system we consider is a single electron in a gate-confined quantum dot (QD) from a 
two-dimensional electron gas (2DEG). The QD is moved along an in-plane trajectory by tuning of the gate voltages. 
The Hamiltonian for the QD-confined electron in the laboratory reference frame is then 

H = H d + H z + Hso + SV(r), (1) 

H d = ^ + u(r-r (t)), ( 2 ) 

Hz = ^g^ B B ■ er, (3) 

HsO = /3-TTyCTx + P+lTxCTy, (4) 

where SV (r) represents a random electric potential, which is always present, whether due to modulation doping or 
barrier disorder. The subscripts d, Z, and SO refer to "dot", "Zecman", and "Spin-Orbit". In Hd, 7r is the electron 
2D momentum (e > 0), given by n = —ifiV + (e/c)A{r) and satisfying [7^, itj] — —ifr^SijkBok', and U (r — Vq (t)) is 
the dot confinement potential with a moving minimum t*o (t) = [xq (t) ,yo (t)]. In this study, we consider a uniform 
linear motion of the QD, where Vq = po/m* = dro {t) jdt is a constant vector. In Hz, Bq is the applied magnetic 
field (with no its unit vector). In Hso- P± = {fi i a), where a and j3 arc the Rashba and Drcssclhaus SO coupling 
constants. The x and y axes are along [110] and [110] directions. If x and y had been defined along [100] and [010] 
directions, the SO term would have taken the usual form Hso = (3(—ir x u x + Ti v Oy) + a(ir x a y — ir y <7 x ). The current 
choice of x and y helps simplify the presentation below. 

The change to the moving reference frame (the rest frame for the electron) is done by a translational transformation 
\tp' (t)) = exp [St (t)] \ip (t)}, where the generator is St (t) = in ■ r (i) /h, and r (t) = [xq (t) ,yo (t)] = (vo x t,vo y t). 
Correspondingly, the Schrodingcr equation becomes 

ih^W{t))=H'W{t)) . (5) 

The new Hamiltonian is 

H' = e ST He- ST + ihd t S T , (6) 

where ihdtST (t) = —tt -po/m*, which we will include in H' d below. If we define A = e St Ae~ ST , then x = x+ [St, x] = 
x + xq (t) and 7T X = tt x + [St, k x ] — tt x — ^Bo z yo (t). After the translational transformation, x and ir x represent the 
position and momentum operators in the moving frame. We re-label them as x' and tt x i, so that x = x' + xq (t) and 
Tt x > = Tt x > — eBo z yo (t) /c. Similarly, we have y = y' + yo (t) and 7ty = ir y > + eBo z xo (t) /c. We can thus write 7r = 
[ir x i — plx', Kyi — pL y >,0], in which pi, (t) = —eBo z /c [—yo (t) , xq (t) , 0] captures the effect of the Lorentz force (note 
that the classical motion of electron under the Lorentz force is given by -# = —e/cvo x B = —e/c-^ (vq x B) ~ jr)- 
The dot Hamiltonian thus takes the form 

H 'd = ^7 + U ( r ') _ {vx'pox + TTy'Poy) /m* 

- { *- p ° f +U(r>)-p L (t).p / m * Pl 



2m* 2m* 



7 



Dropping the c- number terms in H' d , the total Hamiltonian in the moving frame is now 

H' = H' d + H' z + H' so + 5V (r (t) + r') , (7) 

H' d = ^ + P(r'), (8) 

H' z = X - gi x B B-cr, (9) 

H' so = P-n' y ><T x >+P + n' x ><ry, (10) 

where it' = —iliV' + (e/c)A(r') — Pl — Po satisfying [vr^, ir'j] = —ih^SijkBok and B = Bq + AB, where the uniform 
orbital motion induced field is AB = (/3_woy, ft+vox, 0). These are the Eq. (1-4) in the main text. 

CONSTRUCTING THE EFFECTIVE SPIN HAMILTONIAN FOR A MOVING ELECTRON 

As the quantum dot is moved in a semiconductor heterostructure, the spatially random electrical potential 
8V (ro + r') causes the momentum of the QD-confined electron to fluctuate. The electron spin can sense these 
momentum fluctuations via the spin-orbit Hamiltonian (|10|) . and spin decoherence ensues. The QD motion we con- 
sider here is sufficiently slow so that it does not lead to any orbital excitation. This allows us to focus on the effects 
of motion on the spin state of the electron. Below we derive the effective spin Hamiltonian, Eqs. (6) to (8) in the 
main text, by performing a Schrieffer- Wolff transformation and eliminating the orbital degree of freedom. 

The transformation we perform, H = exp(S)H exp(— S), removes the spin-orbit Hamiltonian in the leading order 
by requiring that [H' d + H' Z ,S} = H' so . Formally 



S = 



m— v / 

where the superoperators are defined as h x A = [H' X ,A]. For the harmonic confinement U(r') = ^m*uJ d r' 2 , the 
SO term can be expressed as H' so = ihd {cr ■ £), where £ = (y'/A_, x'/A+, 0) is a vector in the 2DEG plane, in which 
A± = ft,/ (m* f3±) are the spin-orbit lengths. Since [H' d ,H' z ] = 0, the superoperators and l z commute with each 
other. After calculating the commutators, we obtain l z (<r ■ £) = iE z (n x £)-cr, L7 it' = im*r'/h, and L^" 1 [x', y'] = 
—i (hm*ijj d ^) 1 [(tt' x , + m*oj c y') , (n' y , — m*oj c x'^ , where E z = g^sB is the electron Zccman splitting (with ui z = 
E z /h being the Zeeman frequency) and uj c = eBo z / (jn*c) is the cyclotron frequency. Assuming that the Zeeman 
energy is much larger than the SO energy, but much smaller than the orbital excitation energy (m* (/3 2 + a 2 ) <C 
huj z <C hujd), we can expand the generator S only to the first order, 

S = i<T-£ + E z (n-Kl d 1 £)-(T, (12) 
where L^ 1 ^ = Km* J 1 [( n y' — Tn*ui c x'^ /A_, (n' x , + m*ui c y') /A + , 0] . The transformed Hamiltonian is thus 

H" = ihd t S + H' d + H' z + [S, 8V (r)] + • • • , (13) 
in which iftdtS — ^iibABl ■ cr with 

AB L = ^nx — [v Ox /X-,-v 0v /X + ,0]. (14) 
u d g^B 

Recall that AB = (v Oy /\_,v Ox /\ + ,0). Thus AB L = x [\ + AB y /\_, -X_AB X /X + , 0]. As shown in the 

main text, AB < 1 T. Furthermore, uj z <C uJd and w c ~ 1 mcV ~ ujd for a perpendicular 1 T magnetic field in GaAs 
(uj c = for the in-plane magnetic field). Therefore, the correction ABl due to the Lorentz force is generally much 
smaller than AB, and will be dropped from the following considerations. 

The first order term due to the random electric potential is [S, SV (r)] = \g[iBlB x tt(r) ■ cr, in which 

n (r) = [h^t SV (r)] = — ^ [ey (r) /A_, (r) /A+, 0] , (15) 
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where the electric field corresponding to the random potential is e (r) = -'VSV (r) (e > 0). Therefore, the effective 
spin Hamiltonian takes the form 

Heff = \gMB + SB(t))-cr, (16) 
SB(t) = 2BxQ(t), (17) 
□ (t) = (^^(r-oW + OlV), (18) 

where is the orbital wave function. Equations (|T6]) and (|T7| are Eq. (6) and (7) in the main text, which hold in 
the lowest order of spin-orbit interaction and the lowest order of Zeeman splitting (ABl goes to zero in this limit). 
Equation (|17[) shows that at the lowest order, there can be only transverse fluctuations of the effective magnetic field, 
i.e. SB{t) B = 0. 

Equation (fT8|) can be further simplified. Expand the random electric field (r) (i = x', y') around the average QD 
position ro, we obtain 

Si (r) = si (r + r') « Ei (r ) + Ve, (r ) ■ r' + - We* (r ) • r' • r' + ■ ■ ■ . 

Due to the adiabatic condition, the electron always remains in the instantaneous ground orbital state ip(r') = 
exp (— r' 2 /2A 2 ) /X\/tt, up to a magnetic phase that does not affect the calculation of O. Here A~ 2 = 
hr 1 -J '(m*ujd) 2 + (eB z /2c) 2 . For small variation of the gradient of the electric field A 2 We; (ro) <C £-s (f*o) (and 
keep in mind that the average of the linear term with a symmetric ground state wave function vanishes), we retain 
only the zeroth order of £j (r), so that 

n (*) = —^2 (ro (*)) e x> (r (*)) /A+, 0] . (19) 

This is Eq. (8) in the main text. The spatially random electric field that perturbs the orbital motion of the electron 
has now become a temporally random magnetic field that perturbs the rotation of the spin. 



CORRELATION FUNCTION OF THE EFFECTIVE MAGNETIC FIELD 



To calculate the spin relaxation rate of the moving electron, we need to first obtain the temporal correlation 
functions Jij(t) — (SBiSBj (t)) of the random magnetic field that leads to the spin decoherence. Recall that the 
random magnetic field, given by Eq. (jXTJ) and Eq. (|19l) . originates from the spin-orbit interaction and a random 
electric field, the latter from disorder in the substrate material. Thus the temporal correlation of the magnetic 
fluctuations comes from the spatial correlation of the random electric field. Here we choose an isotropic model for the 
random electrical field 

(e< (n) ej (r 2 )> = 6 ij0 if c (Ar/l £ ) , (20) 

where <5y is the Kronecker delta function (i,j = x' or y'), a e is the standard deviation of the electric field, f c (Ar/l e ) 
is the cut-off function as listed in Table I, in which Ar = |n — r%\ is the distance between r\ and and l e is the 
correlation length of the random field. Thus in the moving frame, 

( Si (r (h)) Bj (r (t a ))) = S tJ a 2 f c (\t\ fr c ) , (21) 

where t = ti — t\ and r c = l e /vo is the correlation time. The spatially random electric field in the laboratory frame is 
now a temporally random electric field in the moving frame, which through Eqs. (|17|) and (|19|) becomes a temporally 
random magnetic field. 

The cross product in Eq. (|17[) and the arbitrary direction for the applied magnetic field mean that the magnetic 
correlation is in general quite complex in the (x'y'z) coordinate system we have used so far. To simplify the relaxation 
rate calculations, we first transform to a new XYZ coordinate system, in which we require that (a) Z is along the 
direction of the total magnetic field B and (b) Jij(f) is diagonal in this coordinate system. The first requirement 
dictates that SB is always in the XY plane since SB _L B. This means that 5Bz = and Jzz = 0. The second 
requirement further dictates that the correlation functions are diagonal, so that there are only two independent 
correlation functions Jxx and Jyy- 
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fc(r/le) 


/o(l*l/Tc) 


F c (w) 


1 


exp(-r/Z £ ) 


exp(-|*l /r c ) 


2r c / [1 + w 2 t 2 ] 


2 


1/ [f+r 2 // 2 ] 


1/ [l+t 2 /r c 2 ] 


TTT C exp (— Tc) 


3 


exp (-r 2 /Z 2 ) 


exp (-t 2 /r*) 


y/nT c exp (-oj 2 t 2 /4) 



TABLE I: Fourier transformation of different correlations. 



The XYZ coordinate system can be obtained from the (x'y'z) coordinates by a rotation with Euler angles (p, 9, 
and x- Specifically, first rotate (x'y'z) along the z axis by angle (p to (x"y"z), so that the y" axis is perpendicular to 
the direction of magnetic field n. Then rotate along y" by angle 9 to (x"'y"'Z), so that the Z axis is in the direction 
n. Lastly, rotate along Z axis by angle x- Here angles <p and 9 give the direction of the total magnetic field B in the 
(x'y'z) frame, and x IS determined from the requirement (SBx(0)SBY(t)) = 0. 

After the Euler rotations Rz (x) Ry" (#) R-z (<p)i the field in the XYZ coordinates is given by SB(t) = ~»^ C 

C (t) = [cos x (A x +A v ) + sin X (B x +B y ) , - sin X (A x +A y ) + cos X (B x +B y ) , 0] , (22) 

where A x = — e x > cos (f>/\+, A y = E y > sin^)/A_, 5^ = e x ' cos0sm</>/A + , B y = e y ' cos9coscj)/X-. The condition 
(SBxSByft)) = simply means that (CxCk (*)) = 0. Substituting each component of £ into the equation and 
considering that (ei£j (t)) = 5ij(i^f c (\t\ /r c ), the Euler angle x can be determined as 

2 (A 2 , — A 2 ) n x m y in z , 

tan2 X = ? V + ~ ; / r, (23) 

A 2 , (n 2 , - nlnl, J + A 2 . (n 2 , - n 2 n 2 ,J 

where n is the direction of the magnetic field. 

With the knowledge of all the Euler angles, we can now calculate (CxCx (t)) and ((yCy (t)), 

«x(x(t)) = ([cosx(A x +A y )+sm X (B x +B y )} 2 ) 

(CrCrW) = ([sinx(A x +A y )-cos X (B x +B y )] 2 ) 

where the effective SO length is given by Ajjj 2 = 
2n z / (A+A_). The magnetic correlators are thus 

Jxx(t) 

J YY (t) 

and Jzz(t) = 0, as mentioned earlier. 



= JT°efc (1*1 /Tc) , 

= -La 2 / c (|t|/r c ), 



e (C ± VC 2 - D 2 ) /2, C = (1 - n 2 x ) /A 2 + (f - n 2 ,) /A 2 and D 



2eBa F 



A + m*w^ 
2eBa F 



A_m*u; 2 



/c (1*1/0 



/c (1*1 /Tc) , 



(24) 
(25) 



EFFECTS OF THE DIFFERENT CUT-OFF FUNCTIONS 



The decoherence of the moving electron spin S = cr/2 is governed by the Hamiltonian (|I6[) . In general, the noise 
correlation time r c is much shorter than the spin decay time. In this regime, the dynamics and relaxation of the spin 
is governed by the Bloch equation 0). With purely transverse fluctuations, the longitudinal and transverse relaxation 
rates, 1/Ti and 1/T 2 , are 

^ = |" = JixM + 4 Y (uz), (26) 
J-i ±i 

where the magnetic correlation function in the frequency domain is J^(cj) = g 2 ^f (SBi(0)SBj(t)) cos (ojt) dt. 
Thus 

■#*M = ^^^^ / + °° A (1*1 Ire) cos ( W f) dt, (27) 
(A + m*ujj) J-oo 
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UJTc » 1 


WT C <S 1 


1 


2/ (u; 2 r c ) 


2r c 


2 


TXTc exp (— \bj\ Tc) 


7TT C 


3 


^/tyt c exp (-w 2 r c 2 /4) 





TABLE II: Approximations of F c (oj) in different limits. 



and Jyy(w) is obtained from Eq. (|27[) by substituting A + — > A_. The relaxation rate is then 

2 



1 

Fso 



hoj 2 , 



F so {0,<t>)F c {uj z ), 
) (1 + cos 2 9) + 2a/3 sin 2 9 cos 2q 



(28) 
(29) 



where Fso gives the dependence on the spin-orbit coupling constants and the direction of the magnetic held, and 
F c (lj) is the Fourier transform of f c (\t\ /r c ), as shown in Table I. 




1.0 

Log[v n (m/s)] 



FIG. 4: Spin relaxation rate 1/Ti as a function of the velocity in GaAs QDs for different types of correlation functions in the 
high field limit (in-plane magnetic field B — 1 T, a = and (3 = 300 m/s). 



As shown in Table II, different types of cut-off functions are very similar at the low-field-fast-motion limit, but 
behave dramatically differently in the high-ficld-slow-motion regime. We thus focus on the latter regime, and plot 
the relaxation rate as a function of the QD velocity for different types of correlations in Fig. 01 Overall, 1/Ti is a 
monotonically increasing function of the speed of the QD motion, no matter which type of correlation function is 
used (in particular, 1/Ti is a linear function of vq for the typc-1 correlation function). However, quantitatively the 
relaxation is completely suppressed for the type-2 and -3 cut-off functions because of the exponential suppression from 
exp (— \ljz \ t c ) and exp (— w|t 2 /4) (with a 1 T external field in GaAs and a r c between 1 and 100 ns, we are in the 
limit of tozTc 3> 1). In these cases, spin decoherence is probably dominated by higher-order dephasing processes due 
to the SO coupling 



MEASURING SPIN-ORBIT COUPLING STRENGTH USING CARRIERS TRAPPED BY AN SAW 

A surface acoustic wave (SAW) in GaAs induces a piezo-electric field in the form of E$aw cos (ksAwX — uJsAwt), 
where the k$AW = 2tt/\saw is the wave vector, uisaw is the SAW frequency. The troughs of this propagating 
electric potential can act as a moving QD for electrons, with a velocity at the speed of sound. In the moving frame 
the confinement potential is approximated as Esaw cos (k saw x ) ~ EsawO- ~ k S A W x 2 /2) assuming k$Aw x 1. 
Therefore, the confinement energy can be estimated as u>d = 27ry Eg aw l m tl aw ^ VPs aw / ^saw, where m* = 
0.067m e is the electron effective mass, and Psaw is the RF power that generates the SAW. According to Ref.[12|. 
E S Aw[eV] = 2/25 x iQ(P{dB m ]-23)/20^ We can thus estimate e saw and u d , as shown in Tableim 
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P (dBm) 


Esaw (meV) 


ojo (meV) 


3 


8 


0.5 


13 


25.3 


0.9 


23 


80 


1.6 



TABLE III: Estimation of the confinement energies for different driving power. 



SAW-trapped electrons can help determine the spin-orbit coupling strength in the underlying material. For example, 
in GaAs the electrons and holes that are photo-excited can be picked up by an SAW. The spatial distribution 
of electron spins can then be measured by photoluminescence (specifically polarization of the emitted photons) or 
magneto-optic Kerr rotation [ll, 14). In terms of Kerr rotation, for instance, the spin distribution is expressed as 



9x{d) = 8oe~ d / v ° T2 cos(ujzd/vo) = 6$er d l hB cos(27TK(i), where d is the distance from the origin where the excitons 
are generated, L s = W0T2 is the spin decay length and k = ujz / (27ri>o) is the spatial precession frequency. In the 
absence of an applied magnetic field, the spin precession frequency ujz = g^B^-B /h is completely determined by the 
motion- induced magnetic field AS, 

2m* 

AB = {p_v Qy ,p + v 0x ,Q) . (30) 

This field has been measured to be sizable (25 mT in Ref. [Ill for GaAs), because of the high speed of the QD motion. 
Based on Eq. (|30[) and n = uiz/ {Ittvq) = gfisAB/ (Jivq), one can determine the SO coupling constants by measuring 
the spatial precession frequency k experimentally [6j. This leads to an upper limit of the Dresselhaus SO coupling 
constant at f3 = 300 m/s, which is consistent with other recent experiments However, this /3 value is smaller 

than what was used in earlier theoretical calculations P4ll|. More experimental studies would be needed to clarify 
this issue. 
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